Linking chondrocyte and synovial transcriptional profile to clinical phenotype in osteoarthritis.

OBJECTIVES: To determine how gene expression profiles in osteoarthritis joint tissues relate to patient phenotypes and whether molecular subtypes can be reproducibly captured by a molecular classification algorithm. METHODS: We analysed RNA sequencing data from cartilage and synovium in 113 osteoarthritis patients, applying unsupervised clustering and Multi-Omics Factor Analysis to characterise transcriptional profiles. We tested the association of the molecularly defined patient subgroups with clinical characteristics from electronic health records. RESULTS: We detected two patient subgroups in low-grade cartilage (showing no/minimal degeneration, cartilage normal/softening only), with differences associated with inflammation, extracellular matrix-related and cell adhesion pathways. The high-inflammation subgroup was associated with female sex (OR 4.12, p=0.0024) and prescription of proton pump inhibitors (OR 4.21, p=0.0040). We identified two independent patient subgroupings in osteoarthritis synovium: one related to inflammation and the other to extracellular matrix and cell adhesion processes. A seven-gene classifier including MMP13, APOD, MMP2, MMP1, CYTL1, IL6 and C15orf48 recapitulated the main axis of molecular heterogeneity in low-grade knee osteoarthritis cartilage (correlation ρ=-0.88, p<10-10) and was reproducible in an independent patient cohort (ρ=-0.85, p<10-10). CONCLUSIONS: These data support the reproducible stratification of osteoarthritis patients by molecular subtype and the exploration of new avenues for tailored treatments

A Genome-wide Association Study of Dupuytren Disease Reveals 17 Additional Variants Implicated in Fibrosis

Individuals with Dupuytren disease (DD) are commonly seen by physicians and surgeons across multiple specialties. It is an increasingly common and disabling fibroproliferative disorder of the palmar fascia, which leads to flexion contractures of the digits, and is associated with other tissue-specific fibroses. DD affects between 5% and 25% of people of European descent and is the most common inherited disease of connective tissue. We undertook the largest GWAS to date in individuals with a surgically validated diagnosis of DD from the UK, with replication in British, Dutch, and German individuals. We validated association at all nine previously described signals and discovered 17 additional variants with p ≤ 5 × 10-8. As a proof of principle, we demonstrated correlation of the high-risk genotype at the statistically most strongly associated variant with decreased secretion of the soluble WNT-antagonist SFRP4, in surgical specimen-derived DD myofibroblasts. These results highlight important pathways involved in the pathogenesis of fibrosis, including WNT signaling, extracellular matrix modulation, and inflammation. In addition, many associated loci contain genes that were hitherto unrecognized as playing a role in fibrosis, opening up new avenues of research that may lead to novel treatments for DD and fibrosis more generally. DD represents an ideal human model disease for fibrosis research

Using multivariable Mendelian randomization to estimate the causal effect of bone mineral density on osteoarthritis risk, independent of body mass index

OBJECTIVES: Observational analyses suggest that high bone mineral density (BMD) is a risk factor for osteoarthritis (OA); it is unclear whether this represents a causal effect or shared aetiology and whether these relationships are body mass index (BMI)-independent. We performed bidirectional Mendelian randomization (MR) to uncover the causal pathways between BMD, BMI and OA. METHODS: One-sample (1S)MR estimates were generated by two-stage least-squares regression. Unweighted allele scores instrumented each exposure. Two-sample (2S)MR estimates were generated using inverse-variance weighted random-effects meta-analysis. Multivariable MR (MVMR), including BMD and BMI instruments in the same model, determined the BMI-independent causal pathway from BMD to OA. Latent causal variable (LCV) analysis, using weight-adjusted femoral neck (FN)–BMD and hip/knee OA summary statistics, determined whether genetic correlation explained the causal effect of BMD on OA. RESULTS: 1SMR provided strong evidence for a causal effect of BMD estimated from heel ultrasound (eBMD) on hip and knee OA {odds ratio [OR](hip) = 1.28 [95% confidence interval (CI) = 1.05, 1.57], p = 0.02, OR(knee) = 1.40 [95% CI = 1.20, 1.63], p = 3 × 10(–5), OR per standard deviation [SD] increase}. 2SMR effect sizes were consistent in direction. Results suggested that the causal pathways between eBMD and OA were bidirectional (β(hip) = 1.10 [95% CI = 0.36, 1.84], p = 0.003, β(knee) = 4.16 [95% CI = 2.74, 5.57], p = 8 × 10(–9), β = SD increase per doubling in risk). MVMR identified a BMI-independent causal pathway between eBMD and hip/knee OA. LCV suggested that genetic correlation (i.e. shared genetic aetiology) did not fully explain the causal effects of BMD on hip/knee OA. CONCLUSIONS: These results provide evidence for a BMI-independent causal effect of eBMD on OA. Despite evidence of bidirectional effects, the effect of BMD on OA did not appear to be fully explained by shared genetic aetiology, suggesting a direct action of bone on joint deterioration

Association of a functional microsatellite within intron 1 of the BMP5 gene with susceptibility to osteoarthritis

<p>Abstract</p> <p>Background</p> <p>In a previous study carried out by our group, the genotyping of 36 microsatellite markers from within a narrow interval of chromosome 6p12.3-q13 generated evidence for linkage and for association to female hip osteoarthritis (OA), with the most compelling association found for a marker within intron 1 of the bone morphogenetic protein 5 gene (<it>BMP5</it>). In this study, we aimed to further categorize the association of variants within intron 1 of <it>BMP5 </it>with OA through an expanded genetic association study of the intron and subsequent functional analysis of associated polymorphisms.</p> <p>Methods</p> <p>We genotyped 18 common polymorphisms including 8 microsatellites and 9 single nucleotide polymorphisms (SNPs) and 1 insertion/deletion (INDEL) from within highly conserved regions between human and mouse within intron 1 of <it>BMP5</it>. These markers were then tested for association to OA by a two-stage approach in which the polymorphisms were initially genotyped in a case-control cohort comprising 361 individuals with associated polymorphisms (<it>P </it>≤ 0.05) then genotyped in a second case-control cohort comprising 1185 individuals.</p> <p>Results</p> <p>Two <it>BMP5 </it>intron 1 polymorphisms demonstrated association in the combined case-control cohort of 1546 individuals (765 cases and 781 controls): microsatellite D6S1276 (<it>P </it>= 0.018) and SNP rs921126 (<it>P </it>= 0.013). Functional analyses in osteoblastic, chondrocytic, and adipocytic cell lines indicated that allelic variants of D6S1276 have significant effects on the transcriptional activity of the <it>BMP5 </it>promoter <it>in vitro</it>.</p> <p>Conclusion</p> <p>Variability in gene expression of <it>BMP5 </it>may be an important contributor to OA genetic susceptibility.</p

Combination therapy as a potential risk factor for the development of type 2 diabetes in patients with schizophrenia: the GOMAP study.

BACKGROUND: Schizophrenia (SCZ) is associated with increased risk of type 2 diabetes (T2D). The potential diabetogenic effect of concomitant application of psychotropic treatment classes in patients with SCZ has not yet been evaluated. The overarching goal of the Genetic Overlap between Metabolic and Psychiatric disease (GOMAP) study is to assess the effect of pharmacological, anthropometric, lifestyle and clinical measurements, helping elucidate the mechanisms underlying the aetiology of T2D. METHODS: The GOMAP case-control study (Genetic Overlap between Metabolic and Psychiatric disease) includes hospitalized patients with SCZ, some of whom have T2D. We enrolled 1653 patients with SCZ; 611 with T2D and 1042 patients without T2D. This is the first study of SCZ and T2D comorbidity at this scale in the Greek population. We retrieved detailed information on first- and second-generation antipsychotics (FGA, SGA), antidepressants and mood stabilizers, applied as monotherapy, 2-drug combination, or as 3- or more drug combination. We assessed the effects of psychotropic medication, body mass index, duration of schizophrenia, number of hospitalizations and physical activity on risk of T2D. Using logistic regression, we calculated crude and adjusted odds ratios (OR) to identify associations between demographic factors and the psychiatric medications. RESULTS: Patients with SCZ on a combination of at least three different classes of psychiatric drugs had a higher risk of T2D [OR 1.81 (95% CI 1.22-2.69); p = 0.003] compared to FGA alone therapy, after adjustment for age, BMI, sex, duration of SCZ and number of hospitalizations. We did not find evidence for an association of SGA use or the combination of drugs belonging to two different classes of psychiatric medications with increased risk of T2D [1.27 (0.84-1.93), p = 0.259 and 0.98 (0.71-1.35), p = 0.885, respectively] compared to FGA use. CONCLUSIONS: We find an increased risk of T2D in patients with SCZ who take a combination of at least three different psychotropic medication classes compared to patients whose medication consists only of one or two classes of drugs

Enrichment of low-frequency functional variants revealed by whole-genome sequencing of multiple isolated European populations

The genetic features of isolated populations can boost power in complex-trait association studies, and an in-depth understanding of how their genetic variation has been shaped by their demographic history can help leverage these advantageous characteristics. Here, we perform a comprehensive investigation using 3,059 newly generated low-depth whole-genome sequences from eight European isolates and two matched general populations, together with published data from the 1000 Genomes Project and UK10K. Sequencing data give deeper and richer insights into population demography and genetic characteristics than genotype-chip data, distinguishing related populations more effectively and allowing their functional variants to be studied more fully. We demonstrate relaxation of purifying selection in the isolates, leading to enrichment of rare and low-frequency functional variants, using novel statistics, DVxy and SVxy. We also develop an isolation-index (Isx) that predicts the overall level of such key genetic characteristics and can thus help guide population choice in future complex-trait association studies.Peer reviewe

The Transferability of Lipid-Associated Loci Across African, Asian and European Cohorts

Abstract: Most genome-wide association studies are based on samples of European descent. We assess whether the genetic determinants of blood lipids, a major cardiovascular risk factor, are shared across populations. Genetic correlations for lipids between European-ancestry and Asian cohorts are not significantly different from 1. A genetic risk score based on LDL-cholesterol-associated loci has consistent effects on serum levels in samples from the UK, Uganda and Greece (r = 0.23–0.28, p < 1.9 × 10−14). Overall, there is evidence of reproducibility for ~75% of the major lipid loci from European discovery studies, except triglyceride loci in the Ugandan samples (10% of loci). Individual transferable loci are identified using trans-ethnic colocalization. Ten of fourteen loci not transferable to the Ugandan population have pleiotropic associations with BMI in Europeans; none of the transferable loci do. The non-transferable loci might affect lipids by modifying food intake in environments rich in certain nutrients, which suggests a potential role for gene-environment interactions

Whole-genome sequencing analysis of the cardiometabolic proteome.

The human proteome is a crucial intermediate between complex diseases and their genetic and environmental components, and an important source of drug development targets and biomarkers. Here, we comprehensively assess the genetic architecture of 257 circulating protein biomarkers of cardiometabolic relevance through high-depth (22.5×) whole-genome sequencing (WGS) in 1328 individuals. We discover 131 independent sequence variant associations (P < 7.45 × 10-11) across the allele frequency spectrum, all of which replicate in an independent cohort (n = 1605, 18.4x WGS). We identify for the first time replicating evidence for rare-variant cis-acting protein quantitative trait loci for five genes, involving both coding and noncoding variation. We construct and validate polygenic scores that explain up to 45% of protein level variation. We find causal links between protein levels and disease risk, identifying high-value biomarkers and drug development targets

Accelerating functional gene discovery in osteoarthritis.

Osteoarthritis causes debilitating pain and disability, resulting in a considerable socioeconomic burden, yet no drugs are available that prevent disease onset or progression. Here, we develop, validate and use rapid-throughput imaging techniques to identify abnormal joint phenotypes in randomly selected mutant mice generated by the International Knockout Mouse Consortium. We identify 14 genes with functional involvement in osteoarthritis pathogenesis, including the homeobox gene Pitx1, and functionally characterize 6 candidate human osteoarthritis genes in mouse models. We demonstrate sensitivity of the methods by identifying age-related degenerative joint damage in wild-type mice. Finally, we phenotype previously generated mutant mice with an osteoarthritis-associated polymorphism in the Dio2 gene by CRISPR/Cas9 genome editing and demonstrate a protective role in disease onset with public health implications. We hope this expanding resource of mutant mice will accelerate functional gene discovery in osteoarthritis and offer drug discovery opportunities for this common, incapacitating chronic disease